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Abstract 

Electroweak precision observables, being highly sensitive to loop contributions of new physics, 
provide a powerful tool to test the theory and to discriminate between different models of the 
underlying physics. In that context, the W boson mass, My^r, plays a crucial role. The accuracy of 
the Myy measurement has been significantly improved over the last years, and further improvement 
of the experimental accuracy is expected from future LHC measurements. In order to fully exploit 
the precise experimental determination, an accurate theoretical prediction for Myy in the Standard 
Model (SM) and extensions of it is of central importance. We present the currently most accurate 
prediction for the W boson mass in the Next-to-Minimal Supersymmetric extension of the Standard 
Model (NMSSM), including the full one-loop result and all available higher-order corrections of SM 
and SUSY type. The evaluation of Mw is performed in a flexible framework, which facilitates 
the extension to other models beyond the SM. We show numerical results for the W boson mass 
in the NMSSM, focussing on phenomenologically interesting scenarios, in which the Higgs signal 
can be interpreted as the lightest or second lightest CP-even Higgs boson of the NMSSM. We find 
that, for both Higgs signal interpretations, the NMSSM Myy prediction is well compatible with the 
measurement. We study the SUSY contributions to Myy in detail and investigate in particular the 
genuine NMSSM effects from the Higgs and neutralino sectors. 
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1 Introduction 


Supersymmetry (SUSY) is regarded to be the most appealing extension of the Standard Model (SM), 
as it provides a natural mechanism to explain a light Higgs boson as observed by ATLAS [T] and 
CMS [2]. Supersymmetry realised around the TeV-scale also comes with further desirable features, 
such as a possible dark matter candidate and the unification of gauge couplings. 

The superpotential of the Minimal Supersymmetric extension of the Standard Model (MSSM) 
contains a term bilinear in the two Higgs doublets, 1U(2) ~ ^H 2 Hi. In this term a dimensionful 
parameter, fi, is present, which in the MSSM has no natural connection to the SUSY breaking scale. 
The difficulty to motivate a phenomenologically acceptable value in this context is called the ^-problem 
of the MSSM [3]. This problem is addressed in the Next-to-Minimal Supersymmetric extension of the 
Standard Model (NMSSM), where the Higgs sector of the MSSM gets enlarged by an additional singlet. 
The corresponding term in the superpotential is replaced by a coupling IU( 3 ) ~ XSH 2 H 1 , and the fi 
parameter arises dynamically from the vev of the singlet, S, and may therefore be related to the SUSY 
breaking scale. 

Besides the solution of the //-problem, there are additional motivations to study the NMSSM. The 
physical spectrum contains seven Higgs bosons, which leads to a rich and interesting phenomenology. 
Compared to the MSSM, the singlet field modifies the Higgs mass relations such that the tree-level 
mass of the lightest neutral Higgs boson can be increased. Consequently, the radiative corrections 
needed to shift the mass of the lightest Higgs mass up to 125 GeV can be smaller. This relaxes the 
requirement of heavy stops, or a large splitting in the stop sector, in NMSSM parameter regions where 
the tree-level Higgs mass is larger than the maximal MSSM value (see e.g. Ref. [3]). The NMSSM 
singlet-doublet mixing could also modify the couplings of the 125 GeV boson to explain a potentially 
enhanced rate in the diphoton signal (see e.g. Refs. [SHE]). 

Extensive direct searches for supersymmetric particles are carried out by the LHG experiments. 
These searches have so far not resulted in a signal, which leads to limits on the particle masses, see e.g. 
Refs. [ 9 HII] for a compilation of the results. Indirect methods are complementary to direct searches 
for physics beyond the SM at the LHG and future collider experiments. Whereas direct methods 
attempt to observe traces in the detectors arising from of the direct production of particles of new 
physics models, indirect methods look for the quantum effects induced by virtual exchange of new 
states. Even if not yet seen directly, signs of new physics may show up as small deviations between 
precise measurements and SM predictions. Electroweak precision observables (EWPOs), such as the 
W boson mass, M\y, the sine of the effective leptonic weak mixing angle at the Z boson resonance, 
sin^0^, or the anomalous magnetic moment of the muon, {g — 2)^, (among others) are all highly 
sensitive to loop contributions involving in principle all the particles of the considered model. They 
can both be theoretically predicted and experimentally measured with such a high precision that they 
can be utilised to test the SM, to distinguish between different extensions, and to derive indirect 
constraints on the parameters of a model. Input from indirect methods can be of great interest to 
direct searches for new particles. This was demonstrated, for instance, by the discovery of the top 
quark with a measured mass in remarkable agreement with the indirect prediction [121113] . 

In this paper we focus on the W boson mass. The accuracy of the measurement of Mw has 
been significantly improved in the last years with the results presented by the Tevatron experiments 
CDE [T4] and D0 [T5|. The current world average is [IHlfTT] 

= 80.385 ± 0.015 GeV. (1) 

This precise measurement makes Mw particularly suitable for electroweak precision tests, even more 
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since the precision is expected to be improved further when including the full dataset from the Teva- 
tron and upcoming results from the LHC. Of central importance for the theoretical precision that 
can be achieved on Mw is the top quark mass measurement, since the experimental error on the 
input parameter mt constitutes a dominant source of (parametric) uncertainty, see e.g. Ref. |18] . The 
Tevatron m and LHC |2ni426j measurements of mt have been combined m to yield 

= 173.34 ± 0.27 ± 0.71 GeV. (2) 

In contrast to the Mw measurement, a considerable improvement of the precision on mt beyond Eq. m 
seems less likely at the LHC. Furthermore, it is not straightforward to relate m^^^ measured at a 
hadron collider (using kinematic information about the top decay products) to a theoretically well- 
defined mass parameter. The quantity measured with high precision at the Tevatron and the LHC is 
expected to be close to the top pole mass with an uncertainty at the level of about 1 GeV pHHSO] . 
For the calculation of presented in this paper we adopt the interpretation of the measured value 
mf^^ as the pole mass, but the results could easily be re-expressed in terms of a properly defined 
short distance mass (such as the MS or DR mass). At an e'’"e“ linear collider, the situation would 
improve significantly. Estimates for the ILC show an expected precision ~ 2.5 — 5 MeV and 

Am]^^ = 0.1 GeV [31(132). where the stated precision for mt accounts both for the uncertainty in the 
determination of the actually measured mass parameter and the uncertainty related to the conversion 
into a suitable and theoretically well-defined parameter such as the MS mass. 

For exploiting the precise current and (possible) future Mw measurements, theoretical predictions 
for M\y with comparable accuracy are desired both in the SM and extensions of it. In order to be able 
to discriminate between different models it is necessary that the precision is comparable. In the SM, 
the most advanced evaluation of Mw includes the full one-loop [331134] and two-loop [35H47] result, as 
well as the leading three- and four-loop corrections [48H57j . A simple parametrization for has 

also been developed [58], see also Ref. [59]. Within the SM the LHC Higgs signal at 125.09 GeV [60] 
is interpreted as the SM Higgs boson. Setting ~ 125.09 GeV, the value of Mw can be predicted 
in the SM without any free parameters. The result (for mt = 173.34 GeV, M^ = 125.09 GeV) i^ 

M^ = 80.358 GeV, (3) 

which differs by ~ 1.8 cr from the experimental value given in Eq. (HD- The theoretical uncertainty 
from missing higher-order corrections has been estimated to be around 4 MeV in the SM for a light 
Higgs boson [58] . 

For supersymmetric theories, the one-loop result for Mw |62H73] and leading two-loop corrections 
[MIT] have been obtained for the MSSM. A precise prediction for taking into account 

all relevant higher-order corrections of both SM- and SUSY-type, was presented in Ref. m- A first 
prediction for Mw in the NMSSM has also been presented in Ref. [78]. For the study of other EWPOs 
(mainly focusing on Z decays) in the NMSSM, see Refs. [TSl - lHO] . 

In this work we follow the procedure employed in the MSSM to present a new prediction for 
Mw in the NMSSM with the same level of accuracy as the current best MSSM prediction [61]. We 
combine the complete NMSSM one-loop result with the state-of-the-art SM result and leading SUSY 
higher-order corrections. Our framework allows to output, besides Mw, also the quantity Ar directly 
(see Sect. ED, which summarises all (non-QED) quantum correction to the muon decay amplitude. 
Besides its importance for electroweak precision tests, Ar is needed whenever a theoretical prediction 

^ We updated the SM Mw prediction as discussed below in Sect. m This leads to a small difference compared to 
the SM value given in Ref. m- 
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is parametrized in terms of the Fermi constant (instead of M\y or a{Mz))- Our NMSSM prediction 
for M]y provides the flexibility to analyse SUSY loop contributions analytically and to treat possible 
threshold effects or numerical instabilities. We perform a detailed numerical analysis of Mw in the 
NMSSM with the latest experimental results taken into account. We focus on the effects induced by 
the extended Higgs and neutralino sectors, and in particular on benchmark scenarios where the LHC 
Higgs boson is interpreted as either the lightest or the second lightest CP-even Higgs boson of the 
NMSSM. 

This paper is organised as follows: in Sect. [2]we give a short introduction to the NMSSM, focussing 
on the Higgs and neutralino sectors. In Sect. [3] we describe the determination of the W boson mass in 
the NMSSM. We outline the calculation of the one-loop contributions and the incorporation of higher- 
order contributions. In Sect. 0] we give the numerical results, analysing the NMSSM contributions to 
the W boson mass, before we conclude in Sect. [3 


2 The Next-to-Minimal Supersymmetric Standard Model 


In this section we introduce the NMSSM and specify our notation. We focus on the particle sectors 
which differ from the MSSM. Since the SM fermions and their superpartners appear in the same way 
in both models, the sfermion sector of the NMSSM is unchanged with respect to the MSSM. Also 
the chargino sector is identical to that in the MSSM since no new charged degrees of freedom are 
introduced. For these sectors we use the same notation as employed in Ref. |61j . 

In addition to the two Higgs doublets of the MSSM, the NMSSM also contains a Higgs singlet, S, 
which couples only to the Higgs sectorU Considering the Za-symmetric version of the NMSSM, the 
superpotential takes the form 

VFNMSSM ^ _ dy^QH^ - eyiLHi + XSH 2 H 1 + . (4) 

The new contributions of the Higgs singlet to the soft breaking terms are 

^NMSSM ^ ^MSSM,mod _ ^2 |^|2 _ ^XAxSH 2 Hi + ^ ( 5 ) 

where /;^SSM,mod soft-breaking Lagrangian of the MSSM (see e.g. Eq. (6.3.1) of Ref. [ 8 T]L 

but without the term 6 H 2 FI 1 . The singlet couplings make it possible to dynamically generate an 
effective parameter as 

/Xeff = A(5). ( 6 ) 

The additional contributions (and the modified effective /x term) in the superpotential and in the 
soft breaking terms lead to a Higgs potential which contains the additional soft breaking parameters 
m|, Ax, Af^, as well as the superpotential trilinear couplings A and k. Like in the MSSM, there is no 
CR-violation at tree-level in the couplings of the Higgs doublets. The new doublet-singlet couplings 
allow in principle for CR-violation at tree level, but we will not consider this possibility here. We 
choose all parameters to be real. 

The minimum of the NMSSM Higgs potential triggers electroweak symmetry breaking, after which 
the Higgs doublets can be expanded around their minima according to 


Hi 




(7) 


^For the Higgs doublets and the Higgs singlet we use the same notation for the supermultiplets and for its scalar 
component. 
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Similarly, the singlet scalar component can be expanded as 

S = Vs + ^{4>s + iXs), (8) 

where Vg is the (non-zero) vacuum expectation value of the singlet. 

The bilinear part of the Higgs potential can be written as 

/Xi\ 

Vh = ^ \4>2\+ \ (xi,X2,X5) X 2 1 

K^s) \xsj ( 9 ) 

+ {4>i, 4>2 ) — ’ 

with the mass matrices and M^±^±. 

The mixing of the CP-even, CP-odd and charged Higgs fields occurring in the mass eigenstates is 
described by three unitary matrices , U^, and , where 



These transform the Higgs fields such that the resulting (diagonal) mass matrices become 

KS = {U«y , mS = (f/-’)' and mJ; 5,± = f/^M^*** (t/C)t (n) 

The CP-even mass eigenstates, hi, /12 and h^, are ordered such that rrih^ < ruh^ < ruh^, and similarly 
for the two CP-odd Higgs bosons, oi and 02 - Unchanged from the SM there are also the Goldstone 
bosons, G and G^. Finally, there is the charged Higgs pair, with mass given by 


Ml± = m\ + Ml, - 


( 12 ) 


Here rhA is the effective CP-odd doublet mass given by 

XVs 


m\ = 


{Ax + KVg) . 


(13) 


sin (3 cos /3 

The superpartner of the singlet scalar enlarging the NMSSM Higgs sector is called the singlino, 
S. It extends the neutralino sector with a fifth mass eigenstate. In the basis (P, IT®, iC 2 > <5) the 
neutralino mass matrix at tree level is given by 


M^o = 


(Ml 0 

0 M 2 

—Mz-sw cos P Mzcw cos (3 
sin/3 —Mzcw sin P 

Vo 0 


—Mzsw cos (3 Mzsw sin (3 0 \ 

Mzcw cos 13 —Mzcw sin /3 0 

0 -fieS -Xv2 

-heS 0 -Aul 

-XV 2 -Xvi 2KfieS / 


(14) 


where K = k/X. This mass matrix can be diagonalised by a single unitary matrix N such that 

diag (m^o, m^o , m^o, m^o , m^o ) = N* M^o iVi, 
which gives the mass eigenvalues ordered as m^o < m^o for i < j. 


(15) 
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Figure 1: Left: tree-level diagram with a four-fermion vertex describing muon decay in the Fermi 
model. Right: W boson exchange mediating muon decay in the electroweak SM (in unitary gauge). 

3 Predicting the W boson mass 

3.1 Determination of Mw 

The W boson mass can theoretically be predicted from the muon decay rate. Muons decay to almost 
100% via —7> |82j . This decay was historically described first within the Fermi model (left 

diagram in Fig. [1]). Comparing the muon decay amplitude calculated in the Fermi model to the same 
quantity calculated in the full SM or extensions thereof (the leading-order contribution in unitary 
gauge is depicted in Fig. [1]) yields the relation 

^ + ...,X)). (16) 

This relates the W boson mass to the Fermi constant, G^, which by definition contains the QED 
corrections to the four-fermion contact vertex up to O(a^) |83H87j . and to the other parameters 

Mz and e, which are known experimentally with very high precision. The Fermi constant itself is 
determined with high accuracy from precise measurements of the muon life time |88j . 

The factor Ar in Eq. (jl6l) summarises all higher-order contributions to the muon decay amplitude 
after subtracting the Eermi-model type virtual QED corrections, which are already included in the 
definition of G^. Working in the on-shell renormalization scheme, Eq. (I16p corresponds to a relation 
between the physical masses of the W and Z bosons. 

Neglecting the masses and momenta of the external fermions all loop diagrams can be expressed 
as a term proportional to the Born matrix element |33ll43j 

.'^Loopji — Arj A^Born j ^ ^ Arj . (1^) 

i 

In different models, different particles can contribute as virtual particles in the loop diagrams to 
the muon-decay amplitude. Therefore, the quantity Ar depends on the specific model parameters 
(indicated by the X in Eq. (fT6|) L and Eq. (fTH]) provides a model-dependent prediction for the W boson 
mass. The quantity Ar itself does depend on Mw as well; hence, the value of Mw as the solution of 
Eq. (I16h has to be determined numerically. In practice this is done by iteration. 

In order to exploit the W boson mass for electroweak precision tests a precise theoretical prediction 
for Ar within and beyond the SM is needed. In the next two subsections we describe our calculation 
of Ar in the NMSSM. A new one-loop calculation has been performed which is combined with all 
available higher order corrections of SM- and SUSY-type. 


G, 

V2 
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3.2 One-loop calculation of Ar in the NMSSM 

The one-loop contributions to Ar consist of the W boson self-energy, vertex and box diagrams, and 
the corresponding counter terms (CT). The box diagrams are themselves UV-finite in a renormalizable 
gauge and require no counter terms. Schematically, this can be expressed as 


Ar*-") = W Self-energy -I- W 

( 0 ) 




+ 


Self-energy CT -|- Vertex -|- Vertex CT -|- Box 
-|- Vertex 




J'fw 


(18) 


+ ( 2(5e - 2— + SZu, + t (SZ" + dZ' + SZ"" + SZ"') 


I 


Here S^^(O) denotes the transverse part the W boson self-energy (evaluated at vanishing momentum 
transfer), 5My/ is the renormalization constant for the W boson mass, 5e and 5s^ are the renormal¬ 
ization constants for the electric charge and s* = sin^w, respectively. The 5Z denote different field 
renormalization constants. Since the W boson occurs in the muon decay amplitude only as a virtual 
particle, its field renormalization constant 5Zy/ cancels in the expression for Ar. 

We employ the on-shell renormalization scheme. The one-loop renormalization constants of the 
W and Z boson masses are then given by 


— ^w/z + ^^w/z 1 ^^wiz — ReSj, ^ 


(19) 


where bare quantities are denoted with a zero subscript. The renormalization constant of the electric 
charge is 


eo = (1 -|- Se)e , 6e 


in^^(o) + 


5w s^^(o) 

Cw M| ’ 


( 20 ) 


with 




P 


n^^(o) 




( 21 ) 


Note that the sign ^pearing in front of in Eq. (12011 depends on convention chosen for the SU{2) 
covariant derivativelfl The sine of the weak mixing angle is not an independent parameter in the 
on-shell renormalization scheme. Its renormalization constant 


^w,0 T O'Sw ) 

Sw 


/S?^^(M^) 5]r(M|)\ 

24 V M| ) 


( 22 ) 


is fixed by the renormalization constants of the weak gauge boson masses. 

Finally, the renormalization constant of a (left-handed) lepton held I (neglecting the lepton mass) 
is 

^0 = (l + > (23) 

where denotes the left-handed part of the lepton self energy. 

^We adopt the sign conventions for the SU{2)l covariant derivative used in the code FeynArts | 89H94 | . where (for 
historical reasons) the SU{2)l covariant derivative is defined by — ig 2 l°'WI^ for the SM and 9^ -|- igf2l“lV? 
(N)MSSM. The expressions given here correspond to the (N)MSSM convention. 


7 










Inserting these expressions for the renormalisation constants into Eq. (1181) yields 




+ 2 


S^^(0)-Re[E5r^(M^)] ci; 


rs|^(M|) 




+ (0) - ^ Re 


YWW(-a^2 


(M, 


w> 


Ml 




»W 


^ + Vertex + Box - - (E1(0) + E"(0) + EJ-(O) + E“'(0)) 


The quantity Ar is at one loop level conventionally split into three parts, 


(24) 


= Aa - ^Ap + Ar^ 


(25) 


The shift of the fine structure constant Aa arises from the charge renormalization which contains the 
contributions from light fermions. The quantity Ap contains loop corrections to the p parameter [95], 
which describes the ratio between neutral and charged weak currents, and can be written as 


SP(0) 


(26) 


This quantity is sensitive to the mass splitting between the isospin partners in a doublet [95], which 
leads to a sizeable effect in the SM in particular from the heavy quark doublet. While Aa is a pure SM 
contribution, Ap can get large contributions also from SUSY particles, in particular the superpartners 
of the heavy quarks. All other terms, both of SM and SUSY type, are contained in the remainder 
term Arrem- 

We have performed a diagrammatic one-loop calculation of Ar in the NMSSM according to 
Eq. (1241) . using the Mathematica-based programs FeynArts [89H94] and FormCalc |96] . The NMSSM 
model file for FeynArts, first used in [6], has been adapted from output from SARAH [^1^ . 

The calculation of the SM-type diagrams (being part of the NMSSM contributions) are not dis¬ 
cussed here. This calculation has been discussed in the literature already many years ago msM, and 
we refer to Refs. |43[ I99] for details. The one-loop result for Ar is also known for the MSSM (with 
complex parameters), see Refs. |61[I73|. The calculation in the NMSSM follows along the same lines 
as for the MSSM. However, the result gets modified from differences in the Higgs and the neutralino 
sectors. Below we outline the NMSSM one-loop contributions to Ar, for completeness also including 
the MSSM-type contributions. 

Besides the SM-type contributions with fermions and gauge bosons in the loops (not discussed 
further here), many additional self-energy, vertex and box diagrams appear in the NMSSM with 
sfermions, (SUSY) Higgs bosons, charginos and neutralinos in the loop. Generic examples of gauge- 
boson self-energy diagrams with sfermions are depicted in Fig. [2l Their contribution to Ar is finite 
by itself. The NMSSM Higgs bosons enter only in gauge boson self-energy diagrams, since we have 
neglected the masses of the external fermions. The contributing diagrams are sketched in Fig. [3l These 
contributions are not finite by themselves. Only if one considers all (including SM-type) gauge boson 
and Higgs contributions to the gauge boson self-energy diagrams, the vertex diagrams and vertex 
counterterm diagrams together, the divergences cancel, and one finds a finite result. 

Charginos and neutralinos enter in gauge boson self-energy diagrams (depicted in Fig. |3|), fermion 
self-energy diagrams (depicted in Fig. [5]), vertex diagrams (depicted in Fig. (6] the analogous vertex 
corrections exist also for the other vertex) and box diagrams (depicted in Fig. |7|). The vertex contri¬ 
butions from the chargino/neutralino sector, together with the chargino/neutralino contributions to 
the vertex CT and the gauge boson self-energies are finite. Each box diagram is UV-finite by itself. 












Vi V2 Vi 

fl 


Figure 2: Generic (N)MSSM one-loop gauge boson self-energy diagrams with a sfermion loop; Fl, V 2 = 
7 , Z, and /i, h = i>, I, u, d . 



•Si Si 


Figure 3: Generic NMSSM one-loop gauge boson self-energy diagrams with gauge bosons, Higgs and 
Goldstone bosons in the loop; Fi, F 2 , F 3 = 7 , Z, and si, S 2 = hi, / 12 , h^, ai, 02 , H^, G, G^. 


In order to determine the contribution to Ar from a particular loop diagram, the Born amplitude 
has to be factored out of the one-loop muon decay amplitude, as shown in Eq. (HZl). While most 
loop diagrams directly give a result proportional to the Born amplitude, more complicated spinor 
structures that do not occur in the SM case arise from box diagrams containing neutralinos and 
chargino^. Performing the calculation of the box diagrams in Fig. [7] in FormCalc, the spinor chains 
are returned in the form 

-A^SUSY Box(a) = 

(27) 

-A^SUSY Box(6) = {Uu,<^-Ui,){Uu^U+Ve)b^b) ■ 


The expressions for the coefficients b^a) ^( 6 ) ^.re lengthy and not given here explicitly. In order to 
factor out the Born amplitude 


-FiBorn 


27ra / \ A \ 

2 ..2 [Uel , 


(28) 


the spinor chains in Eq. (I27p have to be transformed into the same structure. We modify the spinor 
chains following the procedure described in Ref. [73] and get 


-FIsUSY Box(a) 
-FI SUSY Box(fe) 


Mb, 


>orn 


(29) 


The coefficients 6 (a) and 6 (f,) contain ratios of mass-squared differences of the involved particles. These 
coefficients can give rise to numerical instabilities in cases of mass degeneracies. In the implemen¬ 
tation of our results (which has been carried out in a Mathematica) special care has been taken of 

"^The same complication occurs in the MSSM and was discussed in Ref. m- 
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Figure 4: Generic NMSSM one-loop gauge boson self-energy diagram with charginos/neutralinos; 
Fi, F2 = 7, = Xt2 and / = X?, 2 ,3,4,5- 





/ N 


l,iy 





Figure 5: Generic NMSSM one-loop fermion self-energy diagram with a chargino/neutralino/sfermion 
contribution; and X° = Xi, 2 ,3,4,5- 


such parameter regions with mass degeneracies or possible threshold effects. By adding appropriate 
expansions a numerically stable evaluation is ensured. 


3.3 Higher-order corrections 


The on-shell renormalization conditions correspond to the definition of the W and Z boson masses 
according to the real part of the complex pole of the propagator, which from two-loop order on is 
the only gauge-invariant way to define the masses of unstable particles (see Ref. |43] and references 
therein). The expansion around the complex pole results in a Breit-Wigner shape with a fixed width 
(fw). Internally we therefore use this definition (fw) of the gauge boson masses. The experimentally 
measured values of the gauge boson masses are obtained using a mass definition in terms of a Breit- 
Wigner shape with a running width (rw). As the last step of our calculation, we therefore transform 
the W boson mass value to the running width definition, to facilitate a direct comparison to the 
experimental value of M\y. The difference between these two definitions is 






2M™’ 


(30) 


where corresponds to the fixed width description, see Ref. [TOO]. For the prediction of the W 
decay width we use 


^ f 

2V27r V 37r y ’ 


parametrized by and including first order QGD corrections. The difference between the fixed- and 
running width definitions amounts to about 27 MeV, which is very relevant in view of the current 
theoretical and experimental precisions. For the Z boson mass, which is used as an input parameter in 
the prediction for M^y, the conversion from the running-width to the fixed-width definition is carried 
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fl 



w 


Figure 6: Generic one-loop vertex correction diagrams in the NMSSM; 
Analoguous diagrams exist for the other W vertex. 



= xt,2 and = Xi, 2 ,3,4,5- 


Figure 7: Generic one-loop box diagrams contributing to the muon decay amplitude in the NMSSM; 
= xt2 and = Xi, 2,3, 4,5- 

out in the hrst step of the calculation. Accordingly, keeping track of the proper definition of the 
gauge boson masses is obviously important in the context of electroweak precision physics. For the 
remainder of this paper we will not use the labels (rw, fw) explicitly; if we do not refer to an internal 
variable, M]y and Mz will always refer to the mass dehnition according to a Breit-Wigner shape with 
a running width (see e.g. Ref. [l3] for further details; see also Refs. |101lll02j l. 

We combine the SM one-loop result (which is part of the NMSSM calculation) with the relevant 
available higher-order corrections of the state-of-the-art prediction for M™. As we will describe 
below in more detail, the higher-order corrections of SM-type are also incorporated in the NMSSM 
calculation of Ar in order to achieve an accurate prediction for For a discussion of the 

incorporation of higher-order contributions to Mw in the MSSM see Refs. [611173] . 

In a first step, we write the NMSSM result for Ar as the sum of the full one-loop and the higher- 
order corrections, 


^^NMSSM ^ ^^NMSSM(o) ^^NMSSM(h.o.) 


(32) 


where /\ 7 -NMSSM(a) NMSSM one-loop contributions from the various particle sectors 


^^NMSSM(q) _ ^^fermion(a) _|_ ^^gauge-boson/Higgs(c 


as discussed in the previous subsection. The term denotes the higher-oder contributions, 

which we split into a SM part and a SUSY part, 


^^NMSSM(h.o.) ^ ^^SM{h.o.) ^^SUSY(h.o.) ^ (-34^ 

The terms describe the SM/SUSY contributions beyond one-loop order. For 

we employ the most up-to-date SM result including all relevant higher-order corrections, while 
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contains the most up-to-date higher-order contributions of SUSY type (see below). The approach fol¬ 
lowed in Eq. ()34p has several advantages. It ensures in particular that the best available SM prediction 
is recovered in the decoupling limit, where all superpartners are heavy, the singlet decouples, and the 
NMSSM Higgs sector becomes SM-like. Furthermore, the approach to combine the most up-to-date 
SM prediction with additional “new physics” contributions (here from supersymmetry) allows one to 
readily compare the MSSM and NMSSM predictions on an equal footing and it provides an appropriate 
framework also for an extension to other scenarios of physics beyond the SM. 

The expression for the higher-order contributions given in Eq. (|34p formally introduces a depen¬ 
dence of the NMSSM result on the SM Higgs mass, which enters from the two-loop electroweak 
corrections onwards. In those contributions we identify the SM Higgs mass with the mass of the 
NMSSM Higgs boson with the largest coupling to gauge bosons. 

For the higher-order corrections in the SM, we incorporate the following contributions 

up to the four-loop order 




(35) 


The contributions in Eq. (I35p consist of the two-loop QCD corrections [35It4nj . the three-loop 

QCD corrections [48H51j . the fermionic electroweak two-loop corrections |42H44j . the 

/ 2 \ 

purely bosonic electroweak two-loop corrections Ar^^^ [45It47j . the mixed QCD and electroweak three- 
loop contributions [5211^ . the purely electroweak three-loop contribution [521153] . 

and finally the four-loop QCD correction [551157] . 

The radiative corrections in the SM beyond one-lom level are numerically significant and lead to 
a large downward shift in M\y by more than 100 MeVo The largest shift (beyond one-loop) is caused 
by the two-loop QCD corrections [35H41] followed by the three-loop QCD corrections [48H51j . 

Most of the higher-order contributions are known analytically (and we include them in this form), 
except for the full electroweak two-loop contributions in the SM which involve numerical integrations 
of the two-loop scalar integrals. These contributions are included in our calculation using a simple 
parametrization formula given in [59] @ This fit formula gives a good approximation to the full result 
for a light SM Higgs (the agreement is better than 0.4 MeV for M\y) [59]. Using this parametrization 
directly for the SM prediction of Ar|“j^ -|- (rather than for the full SM prediction of Mw — an 

approach followed for the MSSM case in Ref. [73] and for the NMSSM case in Ref. [78]) allows us to 
evaluate these contributions at the particular NMSSM value for Mw in each iteration step. The output 
of this formula approximates the full result of Ar^”^ -|- Ar[^“g^ using the fixed-width definition, such 
that it can directly be combined with other terms of our calculation HI In our expression for 
we use the result for given in Ref. ISZj, which contains also contributions from quarks of the 

first two generations and is numerically very close to the result from Ref. [HjH and the result for 


®The corrections beyond one-loop order are in fact crucial for the important result that the Mw prediction in the SM 
favours a light Higgs boson, whereas the one-loop result alone would favour a heavy SM Higgs. 

® In [59] the electroweak two-loop contributions are expressed via AA°‘ ^ -|- = (Aa)^ -I- 2Aa Af^“^ + 

Ar^em\ where a simple fit formula for the remainder term Ar)“m^ is given. The quantity Af^“^ in the second term 
denotes the full one-loop result without the Aa term. 

^It should be noted, however, that the gauge boson masses with running width definition are needed as input for the 
fit formula given in [59]. This is the only part of our calculation where the running width definition is used internally. 

®This is an improvement compared to Ref. |61 |. where the two-loop QCD contributions from Ref. [50] were employed, 
which include only the top und bottom quark contributions. 


12 





Ar("“s) from Ref. [SU]. Both contributions are parametrized in terms of G^. A comparison between 
our evaluation of and the fit formula for given in Ref. [58] can be found in Sect. 14.31 below. 

For the higher-order corrections of SUSY type, see Eq. (|34p . we take the following contributions 
into account, 


^^SUSY(h.o.) ^ ^^SUSY(a2) _ gusY,(aa,) _ SUSY,(a2,aia,(„Q,2) 

rsd q‘2i q‘2i ^ 


(36) 


incorporating all SUSY corrections beyond one-loop order that are known to date. The first term in 
Eq. ([5UD denotes the leading reducible 0{a'^) two-loop corrections. Those contributions are obtained 
by expanding the resummation formula |103j 

1 + Ar =- - -L-^, ( 37 ) 

(1 - Aa) (^1 -h - Arrem 

which correctly takes terms of the type (Aa)^, (A/?)^ and AaAp into accoun1@ if Ap is parametrized by 
Gfj^. The pure SM terms are already included in and because of numerical compensations 

those contributions are small beyond two-loop order m- Thus, we only need to consider the leading 
two-loop terms with SUSY contributions. 


Ar 


SUSY(o2) _ 


red 


= - -^AaAp 


SUSY 


-{Ap' 


SUSYN2 


+ 2:AApSUSY^^; 


.SM 


(38) 


The other two terms in Eq. (|36p denote irreducible two-loop SUSY contributions. The two-loop 
0{aas) SUSY contributions [HIES], contain squark loops with gluon exchange and 

quark/squark loops with gluino exchange (both depicted in Eig. [8|). While the formula for the gluino 
contributions is very lengthy, a compact result exists for the gluon contributions to Ap [74[l75j . Using 
these two-loop results for the SUSY contributions to requires the on-shell (physical) values for 
the squark masses as input. The SU{2) relation implies that one of the stop/sbottom 

masses is not independent, but can be expressed in terms of the other parameters. Therefore, when 
including higher-order contributions, one cannot choose independent renormalization conditions for 
all four (stop and sbottom) masses. Loop corrections to the relation between the squark masses must 
be taken into account in order to be able to insert the proper on-shell values for the squark masses 
into our calculation. This one-loop correction to the relation between the squark masses is relevant 
when inserted into the one-loop SUSY contributions to Mw, while it is of higher order for the two-loop 
SUSY contributions. In our evaluation of Mw this is taken into account by a “mass-shift” correction 
term. For more details see Ref. m- The gluon, gluino and mass-shift corrections, which are identical 
in the MSSM and the NMSSM, are included in our NMSSM result for Mw- 

The third term in Eq. (|36p denotes the dominant Yukawa-enhanced electroweak two-loop correc¬ 
tions to Ap of O(a^), 0{atab) and 0{al) [THllTTl- These contributions consist of heavy quark {t/h) 
loops with Higgs exchange, squark {t/h) loops with Higgs exchange, and mixed quark-squark loops 
with Higgsino exchange, see Fig. [9l The corrections of this kind, which depend on the specific form 
of the Higgs sector, are only known for the MSSM so far [TSIEZ]. It is nevertheless possible to take 
them into account also for the NMSSM in an approximate form. To this end, the considered NMSSM 
parameter point needs to be related to appropriate parameter values of the MSSM. Besides the values 
for tan/3, the sfermion trilinear couplings Aj, and all the soft mass parameters, which can be directly 

®In principle one could also include the term Aa Arseni, which however is numerically small. 
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Figure 8: Generic 0{aas) two-loop self-energy diagrams in the (N)MSSM. Here g denotes a gluon 
and g a gluino; Vi,V 2 = 



Figure 9: Generic 0{a^ + atat + a^) two-loop diagrams, where H denotes a Higgsino and is s either 
an NMSSM Higgs or a Goldstone boson; Vi,V 2 = 


taken over from the parameter point in the NMSSM, we determine the parameters in the following 
way: we set the MSSM /x parameter equal to /Xgff and we use the physical value of the charged Higgs 
mass as calculated in the NMSSM (see below) as input for the calculation of the MSSM Higgs masses. 
This prescription is motivated by the fact that in this way the value of the mass of the charged Higgs 
boson, which is the only Higgs boson appearing without mixing in both models, is the same in the 
NMSSM and the MSSM. The MSSM Higgs masses are calculated with FeynHiggs [lf)4141f)8j . using 
the calculated physical mass value of Mfj± as on-shell input parameter. The MSSM Higgs masses 
and the Higgsino parameter /x determined in this way are then used as input for the calculation of Ap 
to 0{aj), 0{atab), 0{al). In order to avoid double-counting the dominant Yukawa-enhanced elec- 
troweak two-loop corrections in the SM [10911110] have been subtracted according to Eq. (13411 . We find 
that the impact of the Yukawa-enhanced electroweak corrections of SUSY type on Mw is relatively 
small (typically < 1 MeV), and their numerical evaluation is rather time-consuming. In our numerical 
code for Mw in the NMSSM, we therefore leave it as an option to choose whether these contributions 
should be included or not. For the results presented in Sect. 14.41 they have not been included, unless 
stated otherwise. 


4 Numerical results 


4.1 Framework for the numerical analysis 

For the evaluation of the W boson mass prediction, the masses of the NMSSM particles are needed. We 
use the NMSSM on-shell parameters as input to calculate the sfermion, chargino and neutralino masses. 
For the calculation of the Higgs boson masses we use NMSSMTools (version 4.6.0) [rrnmTj f^ For 
other tools that are available to calculate the NMSSM Higgs masses including higher-order radiative 
corrections see Refs. |115H117j . The implementation of the Higgs mass results of Ref. |117| (using 
directly the on-shell parameters as input) is in progress. 


° In two plots below the NMSSM Higgs boson masses at tree-level are used. They are calculated using the tree-level 
relations given in Sect.[2l 
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In NMSSMTools the input parameters are assumed to be DR parameters at the SUSY breaking 
scale. In order to use the NMSSMTools Higgs masses in our result, a transformation from the on- 
shell parameters, needed for our evaluation, to the DR parameters, needed as NMSSMTools input, is 
necessary. This effect is approximately taken into account by transforming the on-shell Xt parameter 
into its DR value by the relation given in Ref. |118] (equations (60) ff in Ref. |118j l. The shift in the 
other parameters is significantly smaller and therefore neglected here. 

We use a setup where the NMSSM parameter space can be tested against a broad set of exper¬ 
imental and theoretical constraints. Besides the constraints already implemented in NMSSMTools^ 
further direct constraints on the Higgs sectors are evaluated using the code HiggsBounds (version 
4.2.0) [120H123] . All programs used for the numerical evaluation are linked through an interface to 
the NMSSM Mathematica code for the W boson mass prediction!^ 

4.2 Theoretical uncertainties 

Before moving on to our numerical results for the W boson mass prediction in the NMSSM, we discuss 
the remaining theoretical uncertainties in the Mw calculation. 

The dominant theoretical uncertainty of the prediction for Mw arises from the parametric uncer¬ 
tainty induced by the experimental error of the top-quark mass. Here one needs to take into account 
both the experimental error of the actual measurement and the systematic uncertainty associated with 
relating the experimentally determined quantity to a theoretically well-defined mass parameter, see 
the discussion above. A total experimental error of 1 GeV on mt causes a parametric uncertainty on 
Mu/ of about 6 MeV, while the parametric uncertainties induced by the current experimental error of 
the hadronic contribution to the shift in the hne-structure constant, Aahadj and by the experimental 
error of Mz amount to about 2 MeV and 2.5 MeV, respectively. The uncertainty of the SM Mw 
prediction caused by the experimental error of the Higgs boson mass, = 0.24 GeV [60], is sig¬ 

nificantly smaller (< 0.2 MeV). In Ref. |32] the impact of improved accuracies of mt and Aohad has 
been discussed. With a precise top mass measurement of Amj = 0.1 GeV (anticipated ILG precision) 
the associated parametric uncertainty in Mw is about 0.6 MeV. 

The uncertainties from unknown higher-order corrections have been estimated to be around 4 MeV 
in the SM for a light Higgs boson {Mjjsm < 300 GeV) [58]. The prediction for M\y in the NMSSM is 
affected by additional theoretical uncertainties from unknown higher-order corrections of SUSY type. 
While in the decoupling limit those additional uncertainties vanish, they can be important if some 
SUSY particles, in particular in the scalar top and bottom sectors, are relatively light. The combined 
theoretical uncertainty from unknown higher-order corrections of SM- and SUSY-type has been esti¬ 
mated (for the MSSM with real parameters) in Refs. [73l[77] as 6Mw ~ (4 — 9) MeV, depending on 
the SUSY mass scale@ Since we include the same SUSY higher-order corrections in our NMSSM cal¬ 
culation as were considered for the uncertainty estimate in the MSSM, the uncertainty from unknown 

NMSSMTools contains a number of theoretical and experimental constraints, e.g. constraints from collider experiments 
(such as LEP mass limits on SUSY particles), B-physics and astrophysics. More details on the constraints included in 
NMSSMTools can be found in Refs. | 1111I119] . 

^^The Mathematica code is linked to a Fortran driver program, calling NMSSMTools and HiggsBounds. The calculation 
of the SUSY particle masses is also included in the Fortran driver. Similarly to the MSSM case m, we plan to 
additionally implement the Mw calculation directly into Fortran in order to increase the speed of the Mw evaluation. 
This will be useful in particular for large scans of the parameter space. 

The lower limit of 4 MeV corresponds to the SM uncertainty, which applies to the decoupling limit of the MSSM. 
For the upper limit of 9 MeV very light SUSY particles were considered. In view of the latest experimental bounds 
from the SUSY searches at the LHC, the (maximal) uncertainty from missing higher orders is expected to be somewhat 
smaller than 9 MeV. 
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ArG) 

ArG“o) 

ArG“«) 




297.17 

36.28 

7.03 

29.14 

-1.60 

1.23 


Table 1: The numerical values (xlO^) of the different contributions to Ar specified in Eq. (j35p are 
given for Mw = 80.385 GeV and = 125.09 GeV. 

higher-order corrections is estimated to be of similar size. 

4.3 SM higher-order corrections 

We compare our evaluation of to the result from the fit formula for given in Ref. |58j . 

In the latest version of Ref. [SH] all the corrections of Eq. (I35p are included. The Mw fit formula 
incorporates the 0{aas) from Ref. [H], whereas we use the 0{aas) from Ref. [37]. These results are 
in good numerical agreement with each other if in both cases the electric charge is parametrized in 
terms of the fine structure constant a. The 0{a‘^as) three-loop corrections included in Eq. (1351) are 
parametrized in terms of We therefore choose to parametrize the 0{aas) contributions also in 
terms of G^. The difference between the G^ parametrization of the QCD two-loop corrections that we 
use here and the a parametrization used in Ref. [58] leads to a prediction for that is ~ 2 MeV 

lower than the result given in Ref. [58]. 

The numerical values of the different SM-type contributions to Ar are given in table [T] for Mw = 
80.385 GeV and Mf^ = 125.09 GeV. The other relevant input parameters that we use are 

mt = 173.34 GeV, = 4.7 GeV, Mz = 91.1876 GeV, Tz = 2.4952 GeV, 

Aaiept = 0.031497686, Aajf^^ = 0.02757, = 137.035999074, 

as{Mz) = 0.1184, G^ = 1.1663787 x lO"*^ GeV^. (39) 

As explained above, the values for the W and Z boson masses given above, which correspond to 
a Breit-Wigner shape with running width, have been transformed internally to the definition of a 
Breit-Wigner shape with fixed width associated with the real part of the complex pole. 

4.4 Results for the Mw prediction in the NMSSM 

We now turn to the discussion of the prediction for Mw in the NMSSM. Our evaluation has been 
carried out for the case of real parameters, consequently for all parameters given in this section the 
phases are set to zero and will not be listed as separate input parameters. 

An earlier result for Mw in the NMSSM was presented in Ref. m- Concerning SUSY two-loop 
contributions, in this result only the part of the contributions to see Eq. (|36h . arising 

from squark loops with gluon exchange is taken into account. As we will show below in the discussion 
of our improved result for Mw in the NMSSM, the two-loop contributions that have been neglected 
in Ref. [78] can have a sizeable impact. A further improvement of our results for the MSSM and the 
NMSSM is that they are based on contributions to Ar that can all be evaluated at the correct input 
value for Mw (using an iterative procedure), i.e. while the evaluation in Ref. [75] makes 

use of the fitting formula for M^ [58]. The corresponding contribution to Ar extracted from the 
fitting formula for M^ is determined at the input value M^ rather than while it is the 

latter that is actually needed for the evaluation in the (N)MSSM (see Ref. [73] for a discussion how to 
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remedy this effect). We have compared our result with the one given in Ref. [7S0 taking into account 
only those contributions in our result that are also contained in the result of Ref. m- We found good 
agreement in this case, at the level of 1-2 MeV on M^. 

Throughout this section, we only display parameter points that are allowed by the LEP limits on 
SUSY particle masses |124j . by all theoretical constraints in NMSSMTools (checking e.g. that the Higgs 
potential has a viable physical minimum and that no Landau pole exists below the GUT scale), and 
have the neutralino as LSP. Unless stated otherwise, we choose the masses of the first and second 
generation squarks and the gluino to be large enough to not be in conflict with the limits from the 
searches for these particles at the LHC0 We make use of the code HiggsBounds [121H123] to check 
each parameter point against the limits from the Higgs searches at LEP, the Tevatron and the LHC. 

4.4.1 Results in the MSSM limit of the NMSSM 

Before turning to the discussion of the genuine NMSSM effects, we show the NMSSM Mw prediction 
in the MSSM limit 


A ^ 0, K ^ 0, K = k/\ = constant, (40) 

with all other parameters (including /ieff) held fixed (such that the MSSM is recovered). In this limit 
one CP-even, one CP-odd Higgs boson (not necessarily the heaviest ones) and one neutralino become 
completely singlet and decouple. In the discussion of in the MSSM limit, the setup for 

the numerical evaluation is introduced and the comparison to the MSSM Mw prediction serves as 
validation of our implementation. 

The left plot of Fig. [10] shows the NMSSM predictions in the MSSM limit (blue curves) as well 
as the MSSM predictions (red curves) for Mw as a function of the stop mixing parameter 
The parameters in Fig. [10] are = 173.34 GeV, tan/3 = 20, /i/eg) = 200 GeV, Mj^^^ = 500 GeV, 
^Qiu /£>12 ~ GeV, Msusy = = M^^ = = 1000 GeV, = A), = At, M 2 = 200 GeV 

and nig = 1500 GeV. For the additional NMSSM parameters we choose ttia = 1000 GeV, A —>• 0, 
K = k/X = 0.5, Ak = —100 GeV (the impact of Af^ on Mw in the MSSM limit is negligible). Here, and 
in the following the prediction for Mw includes all higher-order corrections described above (besides 
the Higgsino two-loop corrections). 

Our approach here is the following: We start from a NMSSM parameter point. We take the 
effective CP-odd donblet mass or the parameter A\ (here rhA = 1000 GeV) as input to calculate 
the NMSSM Higgs boson spectrum. The physical value of the charged Higgs mass (calculated in the 
NMSSM) is used as input for the calculation of the MSSM Higgs masses. As discussed in Sect. 13.31 this 
procedure ensures that the mass of the charged Higgs boson used in onr Mw calculation is the same in 
the NMSSM and the MSSM, since we calculate the MSSM Higgs masses in FeynHiggs (version 2.10.4) 
where the input parameter M^± is interpreted as an on-shell mass parameter. The other parameters 
which occur in both models (tan /3, the sfermion trilinear couplings Af, and the soft mass parameters) 

thank the authors of Ref. [78] for providing us with numerical results from their code. 

The most stringent limits from SUSY searches at the LHC are set on the masses of the first and second generation 
squarks and the gluino, which go beyond ~ 1 TeV. However these limits depend on the model assumptions. Relaxing 
these assumptions, squarks can still be significantly lighter |125| . Substantially weaker limits have been reported for 
the particles of the other sectors, so that third-generation squarks, stops and sbottoms, as well as the uncoloured SUSY 
particles, are significantly less constrained by LHC searches. 

^®The Xt parameter that we plot here is the on-shell parameter. As described in Sect. 14.11 the on-shell value is 
transformed into a DR value, which is used as input for NMSSMTools to calculate the Higgs masses. All numerical values 
given for Xt in this section refer to the on-shell parameters. 
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Figure 10: Comparison of the NMSSM predictions in the MSSM limit (blue curves) for the W boson 
mass (left plot) and the lightest CP-even Higgs mass (right plot) with the MSSM predictions (red 
curves) plotted against the stop mixing parameter Xt. The parameters are given in the text. For 
the two dashed curves (small blue diamonds for the NMSSM predictions in the MSSM limit, and 
red triangles for the MSSM predictions) the tree-level Higgs masses are used. For the solid curves 
(with filled dots) loop-corrected Higgs masses are used: the NMSSM Higgs masses are calculated with 
NMSSMTools, and the MSSM Higgs masses calculated with FeynHiggs. 


are used with the same values as input for the calculation of the physical masses in the MSSM and 
the NMSSM. For the Higgs mass calculation with NMSSMTools the parameter Xt is transformed into a 
DR parameter, while for the calculations its on-shell value is used. The MSSM parameter 

// is identified with the NMSSM effective value 

For the two dashed curves in Fig. [TO] (small blue diamonds for the NMSSM predictions in the 
MSSM limit and red open triangles for the MSSM predictions) the tree-level Higgs masses are used. 
For the solid curve (with filled dots) loop-corrected Higgs masses are used: the NMSSM Higgs masses 
are calculated with NMSSMTools and the MSSM Higgs masses calculated with FeynHiggs. 

The corresponding predictions for the lightest CP-even Higgs mass in the (N)MSSM are displayed 
in the right plot of Fig. fTOl For illustration, in the plots for the Higgs mass predictions the theoretical 
uncertainty on the SUSY Higgs mass is combined with the experimental error into an allowed region for 
the Higgs boson mass, rather than displaying the theoretical uncertainty in the Higgs mass prediction 
as a band around the theory prediction. Consequently, the blue band in the right plot shows the 
region Mjj = 125.09 ± 3.04 GeV, which was obtained by adding a theoretical uncertainty of 3 GeV 
quadratically to the experimental 2cr error. Here Mh represents the corresponding mass parameter 
in the MSSM and the NMSSM (in the considered case Mh in the MSSM and Mh^ in the NMSSM). 
The position of the curves relative to the blue Mh band depends strongly on the other parameters, 
which are fixed here. The range in which the NMSSM parameter points (with NMSSMTools Higgs 
masses) are allowed by HiggsBounds coincides (approximately) with the region in with the lightest 

^^From here on we will leave out the subscript ’eff’ for the /i parameter in the NMSSM 
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Higgs mass is heavy enough to be interpreted as the signal at 125.09 GeV {\Xt\ > 1000 GeV and 
Xt < 1900 GeV). While the tree-level Higgs masses agree exactly in the MSSM and the NMSSM 
in the MSSM limit, we observe a small difference between the masses for the lightest CP-even Higgs 
calculated with FeynHiggs and with NMSSMTools. This discrepancy arises because of differences in the 
higher-order corrections implemented in the two code£l. The tree-level Higgs masses are only used 
in Fig. [To] for illustration. In all following plots (if nothing else is specified) the full loop-corrected 
results for the Higgs masses are used. 

Going back to the left plot of Fig. [TOl we see that the predictions in the MSSM limit 

and the prediction coincide exactly if tree-level Higgs masses are used (which is an important 

check of our implementation). However, using loop-corrected masses, the difference between the 
FeynHiggs and NMSSMTools predictions for the lightest CP-even Higgs mass leads to a difference in 
Mw of ~ 0.8 MeV for small \Xt\. The effect of the difference in the Mw prediction induced by the 
different Higgs mass predictions is contained in the following plots in this section. This should be kept 
in mind when comparing with 

The dependence of the Mw predictions in Fig. [TOjon Xt is influenced both by the loop contributions 
to Ar involving stops and sbottoms, which are identical at the one-loop level in the MSSM and the 
NMSSM, and indirectly via the behaviour of the lightest CP-even Higgs mass. In the chosen example 
the impact of the former contributions is relatively small as a consequence of the relatively high mass 
scale in the stop and sbottom sector. The effect of the higher-order corrections in the Higgs sector is 
clearly visible in Fig. [101 by comparing the full predictions with the ones based on the tree-level Higgs 
masses. As expected from the behaviour of the Mw prediction in the SM on the Higgs boson mass, the 
upward shift in the mass of the lightest CP-even Higgs boson caused by the loop corrections gives rise 
to a sizeable downward shift in the predictions for My/. The local maximum in the M\y predictions at 
about W = 0 is in accordance with the local minimum in the Higgs-mass predictions. The fact that 
the local minima in the Myy predictions are somewhat shifted compared to the local maxima in the 
Higgs-mass predictions is caused by the stop-loop contributions to Ar, whose effect can be directly 
seen for the curves based on the tree-level predictions for the mass of the lightest CP-even Higgs boson 
in the left plot of Fig. [TO] The main contribution of the stop/sbottom sector can be associated with 
Ap and hence depends strongly on the squark mixing. Ap contains terms sensitive to the splitting 
between the squarks of one flavour and terms sensitive to the splitting between stops and sbottoms. 
These two contributions enter with opposite signs, which tend to compensate each other for small and 
moderate values of Xt. 

4.4.2 SUSY higher-oder corrections 

Now we turn to the discussion of the size and parameter dependence of the SUSY two-loop corrections. 
Fig. 111! shows the size of the 0{aas) two-loop corrections. The parameters used here are mt = 
173.34 GeV, tan ,5 = 2, // = 200 GeV, = 1000 GeV, ^ = 1500 GeV, = Ab = 

In NMSSMTools the user can set a flag determining the precision for the Higgs masses. The result from Ref. mi 
containing contributions up to the two-loop level is used if the flag is set equal to 1 or 2, where the two flags correspond 
to the result without (flag 1) and including (flag 2) contributions from non-zero momenta in the one-loop self-energies. 
While in FeynHiggs this momentum dependence is taken into account, we nevertheless And better numerical agreement 
with flag 1 of the NMSSMTools result. For the sake of comparison between the NMSSM and the MSSM predictions 
for Mw it is useful to keep those differences arising from different higher-order corrections in the MSSM limit of the 
Higgs sector as small as possible. We have therefore chosen flag 1 for the Higgs-mass evaluation with NMSSMTools in 
our numerical analyses presented in this paper. As mentioned above, an implementation of our predictions using the 
Higgs-mass evaluation of Ref. is in progress. 
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1000 GeV, M 2 = 600 GeV, rrig = 1500 GeV (solid curves) and nig = 300 GeV (dashed curves), Ax = 
395 GeV, A = 0.57, k = 0.2, = —80 GeV and we vary Msusy = ~ ^ 03 - show the 

results for three values of Xt: Xt = 2 Msusy (left), = 0 (middle) and Xt = —2 Msusy (right). It 
should be stressed here that the parameters for these plots are chosen to demonstrate the possible size 
and the parameter dependence of the SUSY two-loop corrections, however they are partially excluded 
by experimental data: The parameter points in the left plots with Xt = 2 Msusy are HiggsBounds 
allowed for Msusy ^ 800 GeV (apart from a small excluded island around Msusy ~ 550 GeV), 
whereas in the middle and the right plots, the chosen parameters are HiggsBounds excluded for most 
A^susy values. A gluino mass value of nig = 300 GeV is clearly disfavoured by the negative LHC 
search results. Fig. [IT] shows the contribution to the W boson mass, 6 M\y, from the 0{aas) two-loop 
corrections with gluon exchange (dark blue curves), with gluino exchange (orange curves) and from 
the mass-shift correction (pink curves). The shift 5My/ has been obtained by calculating 
twice, once including the corresponding two-loop corrections, and once without, and the two results 
have been subtracted from each other. Starting with the dark blue curves, we find that the gluon 
contributions lead to a maximal shift of ~ 3 MeV in Myr for all three choices of Xt and that the size 
of the gluon contributions decreases with increasing Msusy- Turning to the orange curves, we hnd 
that for nig = 1500 GeV (solid curves) the SMyr shift, induced by the gluino two-loop corrections, is 
small (< 1 MeV) for Vj = 0, while it is up to 3 — 4 MeV for Xt = 2 Msusy and Xt = — 2 Msusy • 
Making the gluino light — choosing nig = 300 GeV (dashed curves) — the gluino corrections can get 
large. For large positive squark mixing, Xt = 2 Msusy, they reach up to 17 MeV for small values of 
A^susy- The gluino corrections can lead to both a positive and a negative Mw shift, depending on 
the stop mixing parameter. Threshold effects occur in the gluino corrections and cause kinks in the 
orange curves, as can be seen in the middle and the right plots. 

The gluon and gluino two-loop contributions are directly related to the mass-shift correction, 
which has to be incorporated in order to arrive at the complete result for the 0 {aas) contributions 
to The pink curves show the impact of this additional correction term. Starting with the 

solid curves {nig = 1500 GeV), we observe that for large stop mixing, Xt = ±2 Msusy, the mass-shift 
corrections are positive and the maximal shift is ~ 4 MeV. For zero mixing the mass-shift corrections 
lead to a large negative shift in Myr (up to —12 MeV for small Msusy)- For nig = 300 GeV, the 
size of the mass-shift correction is smaller. The kinks, caused by threshold effects, can be observed 
(for the same Msusy values) also in the mass-shift corrections. Adding up the gluino and mass-shift 
corrections leads to a smooth curve and no kink is found in the full My/ prediction. This can be seen in 
Fig. [T^ where we plot the sum of the gluon, gluino and mass-shift corrections (all parameters are the 
same as in Fig. [n])- Generally one can see that for large Msusy all contributions decrease, showing 
the expected decoupling behaviour. However contributions from the 0{aas) two-loop corrections up 
to a few MeV are still possible for Msusy = 1000 GeV. 

The Yukawa-enhanced electroweak two-loop corrections of 0{at), 0{atai)), 0{a^) to Ap (“Hig- 
gsino corrections”) in the MSSM can be included in our code, as discussed in Sect. 13.31 To do so, 
we calculate the MSSM Higgs masses as described in Sect. 13.31 (taking the NMSSM charged Higgs 
mass as input for the MSSM Higgs mass calculation) and use them as input for the Ap {0{at), 
0{atai,), 0{a‘^)) formula. The size of these contributions can be seen in Fig. [131 Here, and in 
some of the following plots, we choose modified versions of the benchmark points given in Ref. m, 
which predict one of the CP-even NMSSM Higgs bosons in the mass range of the observed Higgs 
signal, as starting point for our study. Here we take the following parameters: mt = 173.34 GeV, 
tan/3 = 2, n = 200 GeV, M^/^ = 1000 GeV, = 1200 GeV, Mq^ = M^^ = 700 GeV, 

Mp, = 1000 GeV, Ar = Ab = 1000 GeV, Ms = 200 GeV,’"mg = 1500 GeV, aJ = 405 GeV, A = 0.6, 
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Figure 11: Size of the O(aas) two-loop corrections with gluon and gluino exchange. The solid curves 
correspond to nig = 1500 GeV while the dashed curves correspond to nig = 300 GeV. In the left plot 
we set Xt = 2MgusY) in the middle one Xt = 0 and in the right one Xt = — 2MgusY- The plots 
show the contribution to the W boson mass, SM\y, from the O(aas) two-loop corrections with gluon 
exchange (dark blue curves), with gluino exchange (orange curves), and the mass-shift correction (pink 
curves) as a function of MgusY- The parameter points with Xt = 2MgusY are HiggsBounds allowed 
for MgusY ^ 800 GeV, whereas the points with Xt = 0 and with Xt = —2 MgusY predict too low 
Higgs masses and are HiggsBounds excluded for most MgugY values. Note the different scales at the 
y-axis. The parameters used are given in the text. 




Msusy [GeV] WsusY [GeV] 



A4usy [GeV] 


Figure 12: The plots show the full 0{aas) two-loop corrections to Mw (sum of the corrections shown 
separately in Fig. [TT]) as a function of MgusY- The parameters are the same as in Fig. [TTJ The solid 
curves correspond to nig = 1500 GeV while the dashed curves correspond to nig = 300 GeV. In the 
left plot we set Xt = 2MgusY5 in the middle one Vt = 0 and in the right one Xt = — 2MgusY- 


K = 0.18, = —10 GeV, and we vary Xt. These parameter points are HiggsBounds allowed in the re¬ 

gions 700 GeV < Xt < 1000 GeV and 1100 GeV < Xt < 1400 GeV. The left plot shows the NMSSM 
Mw prediction without Higgsino corrections (blue) and including Higgsino corrections (green) plot¬ 
ted against Xt. In the middle plot the shift 5Mw induced by the Higgsino corrections (obtained by 
subtracting the Mw predictions with and withont Higgsino corrections as shown in the left plot) is 
plotted against Xt. We see that the Higgsino corrections can enter the Mw prediction with both 
signs. The numerical effect of the Mw shift, induced by the Higgsino corrections, is relatively small 
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Figure 13: Size of the electroweak ©(af), 0{atah), 0(af) SUSY two-loop corrections. The left plot 
shows the NMSSM Mw prediction without Higgsino corrections (blue) and including Higgsino correc¬ 
tions (green). The middle plot shows the shift 6M-\y induced by the Higgsino corrections (obtained 
by subtracting the M\y predictions with and without Higgsino corrections as shown in the left plot). 
The right plot shows the NMSSM Mw prediction without Higgsino corrections (blue) and including 
Higgsino corrections (green) plotted against the lightest CP-even Higgs mass The black curve 

in the right plot indicates the SM M\y prediction with M^jsm = The grey band indicates the 

1 a region of the experimental W boson mass measurement. The parameters used for these plots are 
given in the text. 


(~ 1 MeV). It was shown in Ref. [77] that the contributions to Mw from the Higgsino corrections can 
be slightly larger (~ 5 MeV) for lighter t/b. The right plot shows the Mw prediction plotted against 
M/jj. We can clearly see here that this scenario, in which the Higgs signal can be interpreted as the 
lightest CP-even NMSSM Higgs, gives a W boson mass prediction in good agreement with the M\y 
measurement indicated by the grey band. 

4.4.3 NMSSM Higgs sector contributions 

Now we turn to the discussion of effects from the NMSSM Higgs sector. In the MSSM the maximal 
value for the tree-level Higgs mass is Mz- One of the features of the NMSSM Higgs sector is 
that the tree-level Higgs mass gets an additional contribution sin^ 2/3, which can shift the 
tree-level Higgs mass upwards compared to its MSSM value (an upward shift can also be caused by 
singlet-doublet mixing, if the singlet state is lighter than the doublet state), and thus reduce the size 
of the radiative corrections needed to ’push’ the lightest Higgs mass up to the experimental value. For 
A = 0.7 and tan /3 = 2 a tree-level value for of 112 GeV is possible [127] . This additional tree-level 
contribution to the Higgs mass, as well as its impact on Mw are shown in Fig. [T4| The parameters 
chosen here are mt = 173.34 GeV, /j = 500 GeV, = 500 GeV, ^q/u/Di 2 ~ GeV, = 

= 1000 GeV, Xt = 2000 GeV, Ar = Ai, = At, Ms = 200 GeV, m-g = 1500 GeV, 
mA = 450 GeV, k = X and = —100 GeV. We vary tan/3 and show the results for different 
values of A. The red curves correspond to the MSSM limit (A —)• 0) while for the other curves the A 
value is given in the corresponding colour. The upper left plot shows the tree-level prediction for the 
lightest CP-even Higgs mass . As expected, the prediction in the MSSM limit approaches its 
maximal value Mz for large tan /3. Increasing A, the prediction decreases for large tan /3, caused 
by doublet-singlet mixing terms. For small tan /3 one clearly sees the positive contribution from the 
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Figure 14: Predictions for Mf^ and Mw as a function of tan/3. The red curves correspond to the 
MSSM limit (A —>• 0) while for the other curves the A values are given in the figure. The upper left plot 
shows the tree-level prediction for the lightest CP-even Higgs mass the upper right plot shows 
Alhi including radiative corrections (calculated with NMSSMTools as described in the text), the lower 
left plot shows the shift 5Mw (calculated as in Eq. (ITTIP from diagrams involving Higgs and gauge 
bosons, and the lower right plot shows the full Mw prediction. The parameters used for these plots 
are given in the text. 


term sin^ 2/3 pushing the tree-level Higgs mass beyond Mz for large A0 The full Mh^ prediction 

The mixing of the h\ state with the heavier singlet leads to a negative contribution to the tree-level Higgs mass, 
which pulls the NMSSM Higgs mass value down (compared to the MSSM case) for intermediate and large tan P values 
(for details see Ref. [4]). At a specific tau/S value this contributions exactly cancels the positive sin^ 2/3 shift at 
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Figure 15: The left plot shows the prediction (blue, solid curve) and the predic¬ 

tion (red) plotted against Xf. In the middle plot, the additional dashed blue curve corresponds to 
_ jVf^(M/jJ -|- M^{Mh) {Mh^ is the mass of the lightest CP-even Higgs of the NMSSM, 
and M/i is the mass of the light CP-even Higgs of the MSSM). The right plots shows the 
prediction plotted against the lightest CP-even Higgs mass The black curve in the right plot 

indicates the SM Mw prediction with Mjjsm = Mf^ ■ The experimental Mw measurement is indicated 
by the grey band; the region Mh = 125.09 ±3.04 GeV is indicated by the blue band. The parameters 
are given in the text. 


(calculated with NMSSMTools as described above) can be seen in the upper right plot. Now we turn to 
the M\\r contributions from the NMSSM Higgs and gauge boson sector, shown in the lower left plot. 
The shift 6M\y displayed here is based on the approximate relation m 


6M\y 


2 



Ar^(“), 


(41) 


where denotes the one-loop contribution from particle sector x (here x=gauge-boson/Higgs), as 

defined for the NMSSM in Eq. (p3]) . The reference Mw value is set here to . The overall 

contribution from the Higgs sector is rather large and negative. As we will discuss in more detail 
below, the Higgs sector contributions here are predominantly SM-type contributions (with M^sm set 
to the corresponding Higgs mass value). The prediction for Mw in the NMSSM is shown in the lower 
right plot. Larger values for M^^ correspond to a lower predicted value for Mw- Thus, for small 
tan/3, where we find a significantly higher predicted value for M^^ for large A than in the MSSM limit 
(arising from the additional tree-level term), we get a lower predicted value for Mw, which is however 
still compatible with the experimental Mw measurement at the 2cr level for the scenario chosen here. 
For tan /? ~ 1 the difference between the W boson mass prediction for A = 0.65 and A —>• 0 is ~ 
25 MeV. The parameter tan/? enters also in the sfermion and in the chargino/neutralino sector. We 
checked that for the parameters used here, the tan /3 dependence of the contributions from these two 
sectors is small compared to the Higgs sector contributions, less than ~ 3 MeV. 

We continue the study of the NMSSM Higgs sector contributions in Fig. [151 In the left plot 
we compare the NMSSM prediction for Mw (blue curve) with the MSSM prediction (red curve). 

the tree level, and the NMSSM Higgs mass value coincides with the MSSM value. In the scenario considered here, this 
happens for all A at the same tan /? value, since we chose k = X. As can be seen in the upper right plot of Fig. 1141 this 
behaviour is approximately retained also in the presence of higher-order corrections in the Higgs sector. 
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The parameters we use here are nit = 173.34 GeV, tan/3 = 2, ^ = 200 GeV, = 1500 GeV, 

^^Q/u/Dt 2 = ^U 3 = ^Qa = 540 GeV, = 1000 GeV, = 1000 GeV, Ms = 

370 GeV, mg = 1500 GeV, A\ = 420 GeV, A = 0.57, k = 0 . 2 , A^^ = —10 GeV, and we vary Xt. The 
NMSSM parameters are allowed by HiggsBounds for Xt > 780 GeV. For Xt > 810 GeV the mass 
of the lightest CP-even Higgs falls in the range of the observed Higgs signal. The MSSM prediction 
is plotted as a comparison to illustrate and discuss the NMSSM effects on Mw- Here (and in the 
following) we do not check any phenomenological constraints for the MSSM parameter point (but 
only for the considered NMSSM scenario). 

The NMSSM prediction for Mw differs from the MSSM prediction by ~ 12 MeV. The chargino/neu- 
tralino contributions can enter with both signs, and we find that in this scenario the relatively small 
/r value causes negative corrections to Ar. On the other hand, small M 2 values tend to give positive 
contributions to Ar. For the chosen parameters, these two effects cancel and contributions from the 
chargino/neutralino sector are very small, 0(0.1 MeV). Gonsequently, different Higgs sector contri¬ 
butions give rise to the difference between the MSSM and the NMSSM curves. Any differences in the 
CP-odd Higgs sector have a negligible impact on the Mw prediction (see also Ref. m)- Since we set 
the charged Higgs masses equal to each other in the two models, differences can only come from the 
CP-even Higgs sector. For this parameter point the second lightest Higgs {M^^ = 150 GeV) has a large 
singlet component (|t7^p ~ 95%), consequently the singlet components of hi and /13 are small. /13 is 
heavy and has no impact on the Mw prediction. Our procedure to calculate the Higgs masses in the 
MSSM and the NMSSM leads to the same charged Higgs masses, but to different predictions for the 
lightest CP-even Higgs masses M^^ and M^. This difference arises from the different relations between 
the charged Higgs mass and the lightest CP-even Higgs mass in the MSSM and the NMSSM. Further 
it also incorporates the (“technical”) difference due to the different radiative corrections included in 
FeynHiggs and NMSSMTools (as analysed above in the MSSM limit). The middle plot of Fig. [T^ shows 
in addition to the NMSSM prediction for Mw (blue) and the MSSM prediction (red), a blue dashed 
curve (with open dots). The dashed blue curve corresponds to — M™(M/iJ -|- M™ (Mi, 1^. 

As one can see the dashed blue curve is very close to the red MSSM curve, thus here the difference 
between the MSSM and the NMSSM Higgs sector contributions to Mw essentially arises from the 
SM-type Higgs sector contributions, in which different Higgs mass values are inserted. It should be 
noted in this context that we have made a choice here by comparing the predictions for a particular 
NMSSM parameter point with an associated MSSM parameter point having the same value of the 
mass of the charged Higgs boson. Accordingly, the predictions for the other Higgs boson masses in 
the two models in general differ from each other, see above, leading to the effect displayed in the left 
plot of Fig. [m Instead, one could have chosen, at least in principle, the associated MSSM parameter 
point such that the masses of the lightest CP-even Higgs masses, M^^ and M^, are equal to each other. 
Also in that case differences in the other parameters in the Higgs sector, including the mass of the 
charged Higgs boson, would induce a shift in the predictions for Mw- 

The right plot of Fig. fTHl shows the prediction plotted against the lightest CP-even Higgs 

mass Mfi^. In this plot we display both the blue band indicating the region M^^ = 125.09 ± 3.04 GeV 
as well as the grey band showing the experimental la band from the W boson mass measurement. 
The black curve in the right plot indicates the SM Mw prediction for M^sm = M^^. It is interesting 
to note that in the NMSSM it is possible to find both the predictions for Mw and for the lightest 
CP-even Higgs mass in the preferred regions indicated by the blue and grey bands in Fig. [151 For the 

^'^The difference in the predictions for the lightest CP-even Higgs masses in the MSSM and the NMSSM, which we 
subtract this way, includes both the difference between the different mass relations in the MSSM and the NMSSM, as 
well as the “technical” difference between the FeynHiggs and the NMSSMTools evaluation. 
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Figure 16: The left plot shows the prediction for (solid curve) and (dashed curve) as a 
function of The region 125.09 ± 3.04 GeV is indicated as a blue band. The middle plot shows the 
singlet components of hi and / 12 , (solid) and (dashed), respectively. The right plot shows the 
prediction. Here the grey band shows the experimental 1 cr band from the W boson mass 
measurement. The parameters used for these plots are given in the text. 


SM, on the other hand, Fig. [TS] shows the well-known result that setting the SM Higgs boson mass to 
the measured experimental value one finds a predicted value for Mw which is somewhat low compared 
to the experimental value. 

Now we want to investigate whether singlet-doublet mixing (a genuine NMSSM feature) has a 
significant impact on the Mw prediction. Such a scenario is analysed in Fig. [161 Our parameters 
are mt = 173.34 GeV, tan/? = 2, ^ = 140 GeV, = 130 GeV ^q/u/Di 2 ~ 1^00 GeV, = 

800 GeV, = 600 GeV, = 1000 GeV, At = 1300 GeV A^ = Ai, = 1000 GeV, M2 = 230 gW, 
rrig = 1500 GeV, Ax = 210 GeV, A = 0.55, k = 0.31, and we vary These parameters are allowed 
by HiggsBounds everywhere apart from —145 GeV < < —105 GeV, and the Higgs signal can 

be interpreted as either hi or /i2. The left plot shows the prediction for (solid curve) and 
(dashed). The corresponding singlet components 17^3 (solid) and (dashed) are shown in the middle 
plot. The third CP-even Higgs is heavy and has a negligible singlet component. For < —120 GeV, 
/i2 is doublet-like and has a mass in the region of the observed Higgs signal (indicated by the blue 
band). In the MSSM, scenarios which allow the interpretation of the Higgs signal as the heavy CP- 
even Higgs involve always a (relatively) light charged Higgs (see e.g. Ref. |128j i. Due to changed 
mass relations between the Higgs bosons, it is possible in the NMSSM to have the second lightest 
CP-even Higgs at 125.09 GeV together with a heavy charged Higgs. Therefore in the NMSSM the 
interpretation of the Higgs signal as the second lightest CP-even Higgs is much less constrained by the 
LHG results from charged Higgs searches [12911130] . The interpretation of the Higgs signal as /i2 in 
this model is always accompanied by a lighter state with reduced couplings to vector bosons. In this 
figure the charged Higgs mass is ~ 280 GeV. For A^ > —100 GeV, hi is doublet-like and has a mass in 
the region of the observed Higgs signal. In the “transition” region (—150 GeV < A^ < —50 GeV) the 
two light CP-even Higgs bosons are close to each other in mass and “share” the singlet component. 
The right plot shows the NMSSM prediction for Mw, which is approximately flat. Accordingly, the 
parameter regions of corresponding to two different interpretations of the Higgs signal within the 
NMSSM lead to very similar predictions for the W boson mass, which are in both cases compatible 
with the experimental result. Even a sizeable doublet-singlet mixing has only a minor effect on the 
Mw prediction in this case. 
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Figure 17: M]y contribution from a light charged Higgs boson. The left plot shows the prediction for 
the CP-even Higgs boson masses in the NMSSM and in the MSSM as a function of the charged Higgs 
mass. The solid curves correspond to the mass of the lightest CP-even Higgs in the NMSSM (blue) 
and the MSSM (red). The dashed curves correspond to the mass of the second lightest CP-even Higgs 
in the NMSSM (blue) and the MSSM (red). The middle plot shows the shift 5M\y (calculated as in 
Eq. (ITTll l induced by the Higgs and gauge boson sector in the NMSSM (blue), in the MSSM (red) 
and in the SM (black) with M^jsm = Mh^. The right plot shows the W boson mass prediction in the 
NMSSM (blue) and the MSSM (red). The parameters used for these plots are given in the text. 


We have demonstrated so far th^ taking Higgs search constraints and the information on the 
discovered Higgs signal into account}^ the genuine NMSSM effects from the extended Higgs sector 
are quite small, and the Higgs sector contributions that we analysed so far were dominated by SM- 
type contributions. This is true in the absence of a light charged Higgs boson, as we will discuss now. 
Light charged Higgs bosons (together with a light CP-even Higgs with small but non-zero couplings 
to vector bosons) can lead to sizeable (non SM-like) Higgs contributions to M\y. This effect can also 
be observed in the MSSM. Although it is not a genuine NMSSM effect, we want to demonstrate the 
impact of such a contribution here. For Fig. [T7]we choose the following parameters mt = 173.34 GeV, 
tan 13 = 9.25, n = 200 GeV, = 300 GeV ^ = 1500 GeV, = 

1100 GeV, At = -2300 GeV = At, = -1500 GeV, M2 = 500 GeV, m-g ="l500 GeV, A = 0.2, 
K = 0.6, Ak = —1370 GeV, and we vary m^. The left plot in Fig. [T71 shows the predictions for the 
masses of the lightest two CP-even Higgs bosons in the NMSSM (blue) and in the MSSM (red) as 
a function of the charged Higgs mass. In both models the second lightest Higgs falls in the mass 
range 125.09 ±3.04 GeV for the chosen parameters. This scenario is essentially excluded by the latest 
charged Higgs searches [12911130] . Nevertheless, we include these plots to illustrate the possible size of 
the contributions from a light charged Higgs. 

The middle plot shows the shift 5Mw calculated as in Eq. (I4T|) with x=gauge-boson/Higgs in the 
NMSSM (blue) and in the MSSM (red) while the right plot shows the full Mw prediction in the 
NMSSM (blue) and in the MSSM (red). As one can see the MSSM and NMSSM contributions to Mw 
are very similar. Since the masses of charginos, neutralinos and sfermions stay constant when varying 
rh/i (or Mjj±), the change in Mw with Mfj± stems purely from the Higgs sector. The Higgs sector 
contribution to Mw comes dominantly from the light charged Higgs, while the lightest CP-even Higgs 

^^Neglecting those experimental bounds one could have very light CP-Higgs bosons with only a small singlet component, 
which would give large contributions to Mw- However this possibility will not be discussed here. 
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gives only a rather small contribution to M\y due to its reduced vector boson couplings. In the middle 
plot the SM result for 5Mw with M^sm = is shown in black. A significant difference between 
the SM Higgs contribution and the MSSM/NMSSM Higgs contributions can be observed. As one can 
see in the right plot, the displayed variation with the charged Higgs boson mass corresponds to about 
a 1(7 shift in M^y. 

4.4.4 Neutralino sector contributions 

We start the discussion of the contributions from the NMSSM neutralino sector, which differs from the 
respective MSSM sector, with Fig. [181 We choose the parameters mt = 173.34 GeV, tan /3 = 3, /i = 
200 GeV, = 1000 GeV, ^ = 1500 GeV, = 650 GeV, = 1000 GeV, 

At = Ar = Ab = 1000 GeV, mg = 1500 GeV, Ax = 580 GeV, A = 0.64, n = 0.25, = —10 GeV, and 

we vary M 2 . In the upper left plot, the blue curve shows the prediction and the red curve the 

prediction. The difference between the NMSSM prediction and the MSSM prediction is small 
for M 2 < 200 GeV and increases for larger M 2 values. The origin of this difference is investigated 
in the other three plots of Fig. ITHl As before our procedure to identify an MSSM point which can 
be compared to the NMSSM point implies different predictions for the lightest CP-even Higgs mass. 
Here we subtract again the difference in the SM contributions, arising from the different Higgs mass 
predictions. The additional blue dashed curve (with open dots) in the upper right plot of Fig. [T8l 
corresponds to — M^(MftJ -|- M^{Mh)- For large M 2 the difference 

between the NMSSM and the MSSM prediction for M\y can be fully explained by the difference 
in the (SM-type) Higgs mass contributions, which arise from inserting different predictions for 
and M}t_. However after subtracting the difference from the Higgs mass contributions we observe a 
sizeable difference between and M^®®^ for small M 2 . This difference stems from different 

sizes of the chargino/neutralino sector contributions between the two SUSY models, which tend to 
compensate the difference between and M^®®^ arising from the Higgs sector. This can be 

seen in the lower left plot, where we display the shift 5My/ (calculated as in Eq. (I4TT) 1 induced by the 
chargino/neutralino contributions in the MSSM (red) and in the NMSSM (blue). At M 2 = 160 GeV 
the chargino mass is 108 GeV and thus just above the LEP limit. The SM\y contribution from the 
chargino/neutralino sector in the MSSM reaches 8.5 MeV in this case@ In the NMSSM the maximal 
6M\y contribution from the chargino/neutralino sector is 16.5 MeV — significantly larger than in the 
MSSM. Both in the MSSM and the NMSSM, the chargino/neutralino contributions decrease when 
increasing M 2 and therewith the chargino and neutralino masses, showing the expected behaviour 
when decoupling the gaugino sector. The largest difference between the NMSSM and the MSSM 
chargino/neutralino contributions is ~ 8 MeV (at M 2 = 160 GeV). The difference arises from the 
neutralino sector, since the chargino sector is unchanged in the NMSSM with respect to the MSSM. 
We will discuss in more detail below why the contributions from the neutralino sector are larger in 
the NMSSM than in the MSSM. The lower right plot of Eig. [18] is similar to the upper right plot, but 
it contains a fourth curve (blue dotted with open diamonds) which was obtained by subtracting the 
different chargino/neutralino contributions, thus it corresponds to 

This curve lies very close to the MSSM prediction. We have therefore identified the contributions 
causing the difference between the and the M^®®^ predictions. 

We continue with the discussion of the neutralino contributions to Myr in the NMSSM in Eig. [TO] 

^^This is not the maximal effect possible for the chargino/neutralino contributions in the MSSM. The 
chargino/neutralino contributions depend on the slepton masses (see diagrams in Figs. I5l7|l . For lighter slepton masses 
the chargino/neutralino contributions in the MSSM can reach up to 20 MeV, as analysed in Ref. |61 |. 
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Figure 18: The upper left plot shows the prediction (blue) and the prediction 

(red) as a function of M 2 . The experimental Mw measurement is indicated as a grey band. The 
upper right plot shows additionally a dashed blue curve (open dots) corresponding to = 

jyjMVisSM _ The lower left plot shows the shift in the W boson mass 6Mw 

(calculated as in Eq. (14111 with x=chargino/neutralino) in the MSSM (red) and in the NMSSM (blue). 
The lower right plot is similar to the upper right plot but it additionally contains the dotted blue curve 
(open diamonds) which corresponds to where 6Mw is the shift 

in Mw induced by the chargino/neutralino contributions. The NMSSM parameter points are allowed 
by HiggsBounds, and falls in the range 125.09 ± 3.04 GeV for M 2 < 725 GeV. The parameters 
used for these plots are given in the text. 


The chosen parameters are mt = 173.34 GeV, tan/3 = 5.5, /r = 200 GeV, Mj^^^ = 245 GeV, 
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Figure 19: The upper left plot shows the masses of the lightest CP-even Higgs boson in the NMSSM 
(blue) and the MSSM (red) as a function of M 2 . The upper right plot shows the prediction for 
(blue) and for (red). The lightest three neutralino masses and the neutralino singlet 

components are displayed in the lower row. The parameters (given in the text) are chosen such that 
the Higgs sectors of the MSSM and the NMSSM are very similar to each other. The parameter region 
in both models is allowed by HiggsBounds and predicts the lightest CP-even Higgs (which is SM-like) 
close to 125.09 GeV. 


^^Q/u/Di2 ^ ^ At = Ar = Ab 1964 GeV, rrig = 

1500 GeV, mA = 1200 GeV, A = 0.62, k = 0.3, A^ = —10 GeV, and M 2 is varied. All parameter 
points are HiggsBounds allowed. Again we get the MSSM prediction by setting the FeynHiggs Mjj± 
input to the value of the charged Higgs mass calculated by NMSSMTools. For this set of parameters 
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Figure 20: The shifts 6Mw in the NMSSM (blue) and in the MSSM (red), calculated taking the 
full chargino/neutralino contribution to Ar into account (solid) and using only the Ap approximation 
(dashed). The parameters are chosen as in Fig. [191 


this procedure leads to a scenario where the MSSM and the NMSSM Higgs boson sectors are very 
similar to each other. Both models predict the lightest CP-even Higgs close to the experimental value 
125.09 GeV, as one can see in the upper left plot of Fig. [19] showing the masses of the two states Mh 
(MSSM, red) and (NMSSM, blue). The difference between and is < 1 GeV, resulting 
in a small (0(1 MeV)) difference in Mw from the Higgs sector contributions. The upper right plot 
of Fig. [T9| displays the W boson mass prediction in the NMSSM (blue) and in the MSSM (red). The 
difference between these two predictions is largest (7 MeV) for M 2 = 150 GeV and (almost) vanishes 
for large M 2 . Since differences in the Higgs sector contributions are quite small, the difference between 
and arises predominately from the differences in the neutralino sector. We note that 

in this scenario both and Mw lie within the preferred regions indicated by the blue and grey 
bands for the whole parameter range displayed in the figure. 

In order to investigate the reasons for the different predictions for the chargino/neutralino con¬ 
tributions we plot the masses of the three lightest neutralino states in the NMSSM (blue) and the 
MSSM (red) in the lower left plot. The other MSSM/NMSSM neutralinos are heavier than 250 GeV 
and hardly affect the Mw prediction. We set here the (unphysical) soft masses Mi and M 2 equal in 
the MSSM and the NMSSM and identify the MSSM p parameter with the effective p of the NMSSM. 
The resulting predictions for the masses of Xi and X 2 are a few GeV lower in the NMSSM than in the 
MSSM. The singlino components of the NMSSM neutralinos, |Vj 5 p, where N was defined in Eq. (|15p . 
are shown in the lower right plot, and we can observe a strong mixing between the five states. The 
singlino components of Xi and X 2 are below 10% for M 2 = 150 GeV and increase up to 40%(20%) for 
Xi(X 2 ) higher M 2 values. The lighter neutralino states (with relatively small singlino component) 
lead to larger contributions from the neutralino sector to M\y in the NMSSM compared to the MSSM. 

In the next step we analyse how well the full Ar contribution of the chargino/neutralino sector can 
be approximated by taking into account only the leading term ^wl^w Ap (defined in Eq. ([26|)). The 
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A/? term contains only the W and Z boson self-energies at zero momentum transfer, thus this approx¬ 
imation neglects in particular the contributions from box, vertex and fermion self-energy diagrams 
containing charginos and neutralinos. The Ap term corresponds to the T parameter of the S, T, U 
parameters |1311I132] , often used to parametrize new physics contribution to electroweak precision ob¬ 
servables. For the plot in Fig. [20] we use the same parameters as in Fig. [T9l Again the blue(red) solid 
curve shows the 5M\y shift as a function of M 2 , calculated as in Eq. (|41l) with x=chargino/neutralino 
in the NMSSM(MSSM) (the two solid curves are identical to the ones in the upper right plot of 
Fig. [T9]) . The two dashed curves show the Mw contributions in the NMSSM (blue) and in the MSSM 
(red) obtained when the full /^ 7 'Chargmo/neutraimo(a) jg approximated by the chargino and neutralino 
contributions to the Ap parameter: 


SMw = -^ 




c2 - 


_A ^charging/neutralino 

SwJ 


(42) 


In the MSSM the Ap term containing charginos and neutralinos provides a very good approximation 
of the full Ar term in the intermediate range 200 GeV < M 2 < 500 GeV. In the range of small and 
large M 2 values, Ap slightly underestimates the full Ar contribution, the difference here is ~ 1.5 MeV 
for M 2 = 150 GeV and ~ 0.5 MeV for M 2 = 750 GeV. In the NMSSM the Ap term gives a 5My/ 
contribution which is larger (> 4 MeV) than the full Ar result for the full M 2 range plotted here. It 
should be noted that the chargino/neutralino sector does not completely decouple for large M 2 in this 
case, which is a consequence of the presence of a light Higgsino, p, = 200 GeV. For M2 = 750 GeV the 
lightest neutralino has a mass of M 2 = 140 GeV, with a singlino component of ~ 40% and a Higgsino 
component of ~ 60%. In this scenario the singlino-higgsino mixing leads to a positive contribution 
to Ap, but to a negative contribution to the Ar terms beyond Ap (we checked that the contribution 
from the box diagrams is negligible for large M 2 values). We also checked that going to large p 
values, the chargino/neutralino sector decouples and all terms vanish. In this scenario the two effects 
largely cancel each other and for large M 2 one finds a small positive value for the full Ar result. 
This however depends on the chosen parameters and the admixture of the light neutralino, e.g. in the 
scenario discussed in Fig. [15] the negative contributions exceed the positive ones so that the full Ar 
result is negative for large M 2 . Thus, we have shown that the Ap approximation for the chargino and 
neutralino contributions works quite well in the MSSM, whereas sizeable corrections to M\y beyond 
the Ap approximation can occur in the NMSSM. 

As a final step we want to discuss the dependence of the Mw prediction in the NMSSM on the 
p parameter, which enters both in the sfermion and in the chargino/neutralino sectors. The left 
plot of Fig. [21] shows the W boson mass prediction in the NMSSM as a function of p, with the 
parameters chosen as mt = 173.34 GeV, tan/3 = 20, Mj^i^ = 250 GeV, ^qiu jbxi ~ 1500 GeV, 
Mq^ = 500 GeV Mjj^ = 1500 GeV, M^,^ = 300 GeV, ^ GeV, At = Ai!'= -2185 GeV, 

AI 2 = 150 GeV, rrig = 1500 GeV, niA = 1500 GeV, A = 0.2, k = 0.6, A^ = —1370 GeV. The 
parameter points are HiggsBounds allowed, and hi falls in the mass range 125.09 ± 3.04 GeV. When 
increasing p, the prediction decreases first, reaches its minimum for p ~ 1100 GeV and then 

rapidly increases. This behaviour can be explained by looking at the contributions to Mw from the 
chargino/neutralino sector (here we take again the full Ar contributions into account) and from the 
stop/sbottom sector. The shift SMw arising from charginos and neutralinos is shown in the middle 
plot of Fig. [21] The chargino/neutralino contribution is largest for small p and decreases with increas¬ 
ing p. Going to larger p the masses of the (higgsino-like) chargino and neutralino states increase and 
the Mw contribution decreases. The shift 6Mw arising from the stop/sbottom sector is shown in the 
right plot of Fig. 1211 The contributions from the stop/sbottom sector (dominated by the Ap contri- 
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Figure 21: Dependence of the W boson mass prediction in the NMSSM on the /r parame¬ 
ter. The left plot shows the prediction, the middle one the 5M\y contribution from the 

chargino/neutralino sector and the right one shows the 5Myi/ contribution from the stop/sbottom 
sector. The parameters are given in the text. 


buttons) get smaller when ^ is increased up to /i ~ 1000 GeV and then start to rise if /i is increased 
further. Increasing //, the splitting between the two sbottoms gets larger (while the stop masses stay 
nearly constant), which implies also an increase of the splitting between stops and sbottoms. The 
counteracting terms in A/? (see the discussion in Sect. 14.4. ip lead to the observed behaviour. 

5 Conclusions 

We have presented the currently most accurate prediction for the W boson mass in the NMSSM, in 
terms of the Z boson mass, the fine-structure constant, the Fermi constant, and model-parameters 
entering via higher-order contributions. This result includes the full one-loop determination and all 
available higher-order corrections of SM and SUSY type. These improved predictions have been 
compared to the state-oUthe-art predictions in the SM and the MSSM within a coherent framework, 
and we have presented numerical results illustrating the similarities and the main differences between 
the predictions of these models. 

Within the SM, interpreting the signal discovered at the LHC as the SM Higgs boson with M^sm = 
125.09 GeV, there is no unknown parameter in the My/ prediction anymore. We have updated the 
SM prediction for Myy making use of the most up to date higher-order contributions. For M^sui = 
125.09 GeV this yields M™ = 80.358 GeV (with a theory uncertainty from unknown higher-order 
corrections of about 4 MeV). The comparison with the current experimental value of = 80.385± 
0.015 GeV shows the well-known feature that the SM prediction lies somewhat below the value that 
is preferred by the measurements from LEP and the Tevatron (at the level of about 1.8 u). The loop 
contributions from supersymmetric particles in general give rise to an upward shift in the prediction 
for Mw as compared to the SM case, which tend to bring the prediction into better agreement with 
the experimental result. 

For the calculation of the Myr prediction, we made use of the highly automated programs FeynArts 
and FormCalc. Our evaluation is based on a framework which was developed in Ref. [6], consisting in 
particular of a NMSSM model file for the program FeynArts and a Fortran driver for the evaluation 
of the masses, mixing angles, etc. needed for the numerical evaluation. The code NMSSMTools is used 
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for the evaluation of the loop-corrected Higgs boson masses. The implementation of another result for 
the NMSSM Higgs masses, obtained in Ref. , is in progress. 

Our improved prediction for the W boson mass in the NMSSM consists of the full one-loop re¬ 
sult, all available higher-order corrections of SM-type, stop-loop and sbottom-loop contributions with 
gluon and gluino exchange of 0{aas), relevant reducible higher-order contributions, as well as an ap¬ 
proximate treatment of the MSSM-type Yukawa-enhanced electroweak two-loop corrections of ©(af), 
0{atab), 0{al). Analytic expressions for all those contributions are implemented, except for the elec¬ 
troweak two-loop contributions of SM-type, for which we make use of the fit formula given in Ref. m- 
The latter allows us to properly evaluate at an NMSSM value for the W boson mass. 

We presented a detailed investigation of the prediction for focussing on the parameter 

regions which are allowed by Higgs searches (tested by HiggsBounds), SUSY searches and further 
theoretical constraints. As a first step we analysed the size of the contributions from stops/sbottoms. 
Since the sfermion sector is unchanged in the NMSSM with respect to the MSSM, we have done this 
study in the MSSM limit, yielding an important check of our NMSSM implementation. We have 
investigated the size of the SUSY two-loop corrections to Mw and found that the 0{aas) corrections 
beyond the pure gluon exchange contributions, which were incorporated in the previous result of 
Ref. [78], can give sizeable contributions. On the other hand, the effect of the Yukawa-enhanced 
electroweak two-loop corrections of 0{a‘^), 0{atab), 0{al) stays numerically relatively small in the 
allowed parameter region. 

Concerning the investigation of genuine NMSSM effects, we started our discussion with the Higgs 
sector contributions to Mw. The tree-level prediction for the lightest CV-even Higgs mass is modified 
by an additional term in the NMSSM as compared to the MSSM, which (for small tan /3) leads to an 
upward shift of the tree-level Higgs mass. Therefore, in that region, the radiative corrections needed 
to push the Higgs mass to about 125 GeV can be smaller than in the MSSM, which implies that 
lighter stop masses and a smaller stop mixing are possible. We investigated a scenario where this 
additional tree-level term gives rise to a higher prediction than in the MSSM limit. The impact 
on the Mw prediction is a downward shift (of ~ 25 MeV in the considered example) as compared 
to the corresponding prediction in the MSSM. In the prediction for Mw this contribution from the 
Higgs sector enters together with other SUSY loop contribution to Ar yielding an upward shift in 
Mw compared to the SM. The overall effect is such that also in a scenario of this kind a very good 
agreement between the theoretical prediction and the experimental result can be reached. We have 
furthermore investigated the effect of doublet-singlet mixing. While a sizeable doublet-singlet mixing 
can occur in the region where the two NMSSM Higgs states hi and /12 are close to each other in mass, 
we find that it has only a minor effect on the Mw prediction. 

In the NMSSM the Higgs signal seen at the LHC can be interpreted both as the lightest and the 
second-lightest CV-even Higgs boson of the spectrum. Both interpretations give predictions for the 
W boson mass in good agreement with the M\y measurement. In the NMSSM the interpretation of 
the LHC signal as the second-lightest CV-even Higgs /12 is possible together with a relatively heavy 
charged Higgs. This is different from the situation in the MSSM, where all Higgs states have to be 
light in this case, so that such a scenario can be probed by searches for charged Higgs bosons in top- 
quark decays. As a consequence, the interpretation of the observed Higgs signal as the second-lightest 
CV-even Higgs boson is much less constrained in the NMSSM compared to the MSSM. 

For completeness, we have nevertheless briefly investigated also the case of a light charged Higgs 
boson. We have found that a light charged Higgs boson (together with a light CV-even Higgs with 
reduced but non-zero couplings to gauge bosons) can in principle give very signihcant contributions 
to M\y (as in the MSSM). In that case large deviations from the SM Higgs sector contributions occur. 
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but as discussed above scenarios of this kind are severely constrained by limits from charged Higgs 
searches at the LHC. Generally we find that taking all available constrains on the Higgs sector into 
account, the specific NMSSM effects of the Higgs sector to Mw are relatively small. 

On the other hand, the extended neutralino sector of the NMSSM can lead to a sizeable difference 
between the W boson mass predictions in the NMSSM and the MSSM. The chargino/neutralino 
contributions to My/ can be larger in the NMSSM compared to the MSSM, where in the scenario 
which we studied the difference reaches ~ 8 MeV. Assuming the same values for the soft mass 
parameters in the MSSM and the NMSSM and choosing /r = /Ueff) the mixing with the singlino leads 
to shifts in the neutralino masses as compared to the MSSM case. In the considered scenarios the 
lightest NMSSM states turned out to be lighter than the corresponding MSSM states. They also have 
a relatively small singlino component, which causes the resulting contributions to the prediction for 
Myy to be larger than in the MSSM. While light wino/bino states typically give positive contributions, 
light higgsinos can give contributions entering with both signs. 

As a final step of our analysis, we compared the My/ prediction calculated with the full Ar 
to the one where the full result is approximated by the contribution to Ap. We found that the 
difference between the full result and the Ap approximation can be sizeable in the NMSSM, where 
the approximation can lead both to an over- or an underestimate of the full result. Light neutralinos 
with a significant higgsino-singlino mixing tend to give a positive contribution to Ap, but a negative 
contribution to the Ar terms beyond Ap. It therefore depends on the exact patterns of the admixture 
with the singlino whether the neutralino sector of the NMSSM leads to an upward or downward shift 
in the prediction for M/y in comparison with the MSSM. 

We have demonstrated that the prediction for the W boson mass arising from the relation with 
the Z boson mass, the Fermi constant and the fine structure constant in comparison with high- 
precision measurements of those quantities provides a high sensitivity for discriminating between the 
SM and possible extensions of it. With further improvements of the experimental accuracy of Myy, 
possible improvements in the determination of and further information on possible mass spectra 
of supersymmetric particles - either via improved limits or the discovery of new states - the impact 
of this important tool can be expected to be even more pronounced in the future. 
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